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I. INTRODUCTION 

Complex networks have been attracting great attention 
for decades. They may describe many real-world systems 
in the social sciences, biology, computer science, telecom- 
munication and others [HE|. Mathematical description 
of networks is provided by the graph theory Q . Graph is 
a set of vertexes (nodes) connected by edges (links). The 
main local characteristic of a graph is the node degree, i.e. 
the number of links incoming to or outgoing from a node. 
For almost fifty years, the paradigm of "typical node" has 
been present in the science of networks. Networks of typ- 
ical nodes were described by Erdos and Renyi Q (classi- 
cal random graphs — CRG). In their model, iV nodes are 
connected randomly with L links: each inter- node link is 
realized with the probability p = 2L/[N(N — 1)]. In this 
model, the node degree distribution is given by a Pois- 
son law, i.e. Pk{k) = exp(— {k}){k} k /k\, where {■ ■ • } 
denotes the mean over all N nodes, and the node degrees 
observed on a CRG fluctuate around {k}. 

As pointed out by Albert and Barabasi in their semi- 
nal paper networks in real world more often exhibit 
a power-like degree distribution, i.e. Pk(k) oc k^ 1 . In 
the Albert-Barabasi (A-B) model, the node degrees as- 
sume all integer values in the thermodynamic limit and 
there is no characteristic value of the degree. Thus, with 
this observation the Hungarian mathematicians' world of 
networks with typical nodes became a world of scale-free 
networks. 

CRG and A-B networks are two examples belonging to 
two different families of networks 6j. The first one be- 
longs to the so-called homogeneous networks, which may 
be described via a statistical ensemble. The A-B net- 
works have temporal structure, as they come into being 
via growth process. The A-B network is an example of 
causal network. 
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For networks, the act of growing means subsequent at- 
tachments of new nodes, each with M links, to previously 
existing nodes. The procedure of selection of those "old" 
nodes influences the network topology and the degree dis- 
tribution. When old nodes are selected randomly — i.e. 
the probability of attachment is the same for all nodes 
— exponential networks are created and the nodes degree 
distribution is an exponential one 0] . On the other hand, 
when the attachment is preferential — i.e. the probabil- 
ity of choosing a node is proportional to its degree - 
the degree distribution is power-like and network can be 
termed as scale-free 

The number of edges M also influences the network 
topology: 

• when M = 1, the path between any pair of nodes 
is unique; the growing structure is called a tree, 

• when M > 1, cyclic paths are possible and graph 
looses its tree-like properties, 

• when M > 1 and chosen old nodes all are different, 
multiple edges in the network are absent and the 
structure is a simple graph. 

Such attaching procedure prevents possibility of loops, 
i.e. self-links. 

Several characteristic of real or simulated networks 
may be practically useful. For example many papers 
discuss the networks resistance to possible damage H, 
their tolerance on random and/or intentional attack [8| 
or transport properties in terms either of the percolation 
theory j^] or of the shortest path finding 0, ^| . New- 
man et al. applied the generating function formalism 12] 
to evaluate the number of nodes 

~ _ 2 — m m — 1 (i\ 

in subsequent (m-th) layer from a randomly chosen ori- 
gin 0| . In Eq. z\ and z^ are typical values of the 
number of nodes nearest neighbors and the number of 
nodes next-nearest neighbors, respectively. The first one 
(zi) is obviously equal to average node degree z\ = {fci}. 
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The latter (22) was evaluated lately by Shargel et al. 01 



as 



Z2 = {k 2 } - {k}. 



(2) 



Basing on the same technique, Motter et al. [l5| de- 
rived the average node-to-node distance dij dependence 
on the product kikj of the node degrees for random net- 
works: 



(ckj) =A-Bhx(kikj), 



(3) 



where (■ • • } denotes the average over all node pairs, the 
product of the pair degrees being equal to kikj . Recently, 
Holyst et al. |l6j have confirmed this dependence numeri- 
cally and presented some examples of real- world networks 
which obey Motter et al. theoretical predictions. 

In this paper we check if Motter et al., Holyst et al. 
and Shargel et al. predictions apply to the growing expo- 
nential networks. Namely, we evaluate number of neigh- 
bors in subsequent layers. The node-to-node distance vs. 
product of their degrees is also simulated. For complete- 
ness, the calculations and discussion include the growing 
scale-free A-B networks. 

In the next section we explain our numerical approach. 
In section IIHI we present results of Monte Carlo simu- 
lations of the average number of subsequent neighbors 
(|III All and the inter-nodes distance dependence on the 
product of their degrees IjHI 
marizes the results. 



The last section sum- 



II. NUMERICAL APPROACH 

Numerical approach is based on an "on-line" construc- 
tion of the distance matrix D during the network growth 
n element dij of the distance matrix 
gives the length of the shortest path between nodes i and 
j, i.e. the minimal number of edges which connect these 
vertexes. The numbers rfy in i-th row/column inform 
how far is the node i from another node j. Then, the 
number z m (i) representing the number of occurrences of 
the m value in the i-th matrix row/column gives the in- 
formation how many neighbors of the node i are at the 
distance m po| . The average number of the matrix el- 
ements of a given value in all rows/columns — i.e. in 
the whole matrix — gives a typical number z m of subse- 
quent neighbors, for example, the nearest neighbors for 
m = 1, the next-nearest neighbors for m = 2, the next- 
next-nearest neighbors for m = 3, etc. Additionally, the 
number of unities ("1") in the i-th row/column gives de- 
gree of the i-th node: k(i) = z\(i). 



III. RESULTS OF SIMULATIONS 

We construct the distance matrix D for N = 
nodes. The results are averaged over N Iun = 10 4 
dependent simulations. 
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A. Number of nodes in subsequent layers 
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Fig. shows how the deviation S m = z\ "' z' 2 
z m between z m calculated from Eq. (JIJ and from the 
direct simulation behaves as the function of the system 
size N for m — 3, 4, 5. As one can see, starting with N w 
100 this difference decreases with N for the exponential 
networks. For the scale-free networks either the number 
of nodes (N = 10 3 ) is still too small to observe a good 
agreement between z m and z 2 ~ m z 2 n ~ 1 or Eq. (JIJ docs 
not hold for the A-B graphs. 

By construction, the average number of the nearest 
neighbors z\ is 2 and 4 for trees and simple graphs, re- 
spectively. The number of next-nearest neighbors z-i de- 
pends on the applied rules of growth: when the growth 
is governed by the preferential attachment rule, we have 
Z2 ~ 14 and Z2 ~ 38 for M = 1 and M = 2, respectively. 
For the exponential networks, these numbers are Z2 ~ 4 
(M = 1) and Z2 » 17 (M = 2). As the average number 
of the nearest neighbors z\ is exactly equal to the average 
nodes degree {k}, it may be evaluated from the degree 
distribution Pk(k) as {k} — kPk(k), as well. For 

the exponential network this distribution [3, H(j is given 
by 



for M = 1, 



3/4-(3/2)- fe forM = 2, 



Pk(k > M) = 

while for the scale-free networks |2jj, |22j it is 

2M(M + 1) 



Pk(k>M) = 



(fc + 2)(fc + !)£;■ 



(4) 



(5) 



The mean number of the next-nearest neighbors (22) 
may be evaluated basing on Eq. which diverges for 
scale- free networks with N — > 00. For finite but large 
network this sum 



N-l 



N-l 



a = k(k-l)P k (k) = 2M(M+1) ^ 



1 



k=M 



k=M 



(k + 2)(k + 1) 



(6) 

grows logarithmically with N, a — 3.99 In A — 7.55 (M = 
1) and a = 11.96 In AT— 22.6 (M = 2) as presented in Fig. 
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The results are collected in Tab. [I] 



B. Node-to-node distances and node degrees 

Using t he g enerating function formalism 0, Mot- 
ter et al. [15| derived an expression for the length of the 
shortest path between the nodes for a given value of the 
product of connectivities kikj: 



(dij) = 



1 + 



\n(N Zl ) 
\n(z 2 /zi) 





1 




\n(z 2 /zi)_ 



ln(fcjfcj) = 



A Xh - B Xh \u(kikj). 



(7) 



3 



scale-free tree (M=1) 



exponential tree (M=1) 




FIG. 1: Dependence 8 m = z 1 2 - m z 2 m - 1 - z m (m = 3,4,5) on network size N for growing exponential and scale-free trees 
(M = 1) and simple graphs (M = 2). 



TABLE I: Average number of the nearest neighbors z\ and 
the next-nearest neighbors z 2 for different evolving scale- 
free and exponential networks with N = 10 3 nodes. The 
results are averaged over N Tun — 10 4 samples. Theoreti- 
cal predictions of the average nodes degrees {&}(= zi) and 
{k 2 } — {k}(— z 2 ) are also included. Four last lines show 
the least-square fit coefficients A and B in the dependence 
(dij) = A — B\n(kikj) and their predictions A th and B th 
given by Eq. Q. 
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FIG. 2: Dependence a(N) for growing scale-free networks. 



Lately, such a kind of dependence (dij) vs (kikj) has 
been shown to be valid in few real-world networks, in- 
cluding biological and scientific papers citation networks, 
public-transportation systems of several Polish towns, 
and simulated CRG and A-B networks [lfij . 

Here we show that this logarithmic dependence 
holds for exponential networks with M = 1,2. The re- 
sults are presented in Fig. |3| and in Tab. [I] The least- 
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square fit was confined to the two first decades of kikj 
values for the scale-free networks and to the one-and-half 
decade of kikj values for the exponential ones. 

IV. DISCUSSION AND CONCLUSIONS 

Generating function mechanism 0, 0, 0> EH Eil has 
a mean-field nature and should work only for homoge- 
neous trees. This mechanism is founded upon the as- 
sumption that there are no correlations between nodes 
degrees. But this assumption does not hold for growing 
(causal) trees, where the oldest nodes — probably well 
connected — are geometrically close |27j] . 

However, the Motter et al. formula JJJ works surpris- 
ingly nice also for growing networks, where triangles and 
other cyclic paths are possible. The Motter et al. predic- 
tions agree with simulation particularly well for exponen- 
tial networks, where also z m were reproduced quite fairly. 
For scale- free networks this agreement is only qualitative. 
It seems, that theoretical predictions given by Eqs. 
and 10 (and obviously given by Eq. agree with re- 
sults of the simulations for the networks for which degree 
distribution gives finite {k 2 }. 

Average number of vertexes in all generation z m is well 
known for homogeneous [2{| and causal [24| trees. The 



number z m of m-th neighbors derived in |l3l Il5| agree 
very well for small m = 3,4,5 with the results of com- 
puter simulation for exponential networks where the old 
nodes, to which the new nodes are being attached, are 
chosen randomly. 

Again, this should be valid for trees, but it works also 
nicely for M = 2 when random attachment is used. On 
the other hand the sum YlkLM^ik — l)-Pfc(^) diverges 
for power-like distributions Pfc(fc). For finite but large 
lattices this sum (a, Eq. ©) increases logarithmically 
with the system size A. In all four investigated cases 
simulated number Z2 (given by number of "2" in distance 
matrix) agrees with {k 2 } — {k} (averaged over all graph 
nodes). 

For larger m formula fails when applied to real 
networks, i.e. with finite A. Usually, the second layer 
contains more nodes than the first one, from which fol- 
lows that Z2 > Z\ . Then — accordingly to Eq. (0 — z m 
increases with m G Z, but for finite systems it must start 
to decrease for large m. In particular, any of A nodes 
which constitute the network has no neighbors in A-th 
layer and does not posses any A-th neighbors (z m = for 
m > A) . The distribution of the node-to-node distances 
for the growing networks discussed here were shown in 
R^, |2(j and evaluated analytically in case of trees in Ref. 
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Still, the method of evaluation of z m (m = 3,4,5) 
based on Eq. (|TJ may be quite useful. The main effort 
should be paid to a theoretical evaluation of the average 
number of nodes in the second layer, i.e. the number of 
occurrences of "2" in the distance matrix, basing only on 
the degree distribution Pk{k). Such an evaluation of Z2 
would allow, in principle, to reproduce the whole function 

Although the node-to-node distance (dij ) depends log- 
arithmically on the product of the node i and j degrees 
(Fig. | Eq. ©, Ref. 0), the dependence of the to- 
node distance on the node degree is not a trivial one . 
We have demonstrated, that Eq. (J3J can be extended to 



the case of the growing exponential networks. The Mot- 
ter et al. predictions of values A and B in Eq. J3J given 
by Eq. JJJ agree for these networks quite fairly. 
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